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Abstract 

Air shower experiments have detected cosmic ray events of energies upto 300 EeV. 
Most likely these cosmic rays have originated from compact objects. Their exact sources 
are yet to be identified. It has been suggested before that gamma ray bursts are 
possible sources of ultra- high energy cosmic rays. The two models of gamma ray burst 
emissions most often discussed are the internal and external shock models. We have 
calculated the proton spectrum above 60EeV from all gamma ray bursts distributed 
upto a redshift of 0.02 in the internal shock model assuming redshift and luminosity 
distributions consistent with observations, log normal distributions for their values of 
Lorentz factors, variability times and duration of bursts. Within the external shock 
model we have calculated the proton flux above 60EeV from all nearby gamma ray 
bursts assuming log normal distributions in their values of total energies, Lorentz 
factors at the deceleration epoch and compared with the observed data. We find that 
gamma ray bursts can produce cosmic ray proton flux comparable to the flux observed 
by the Pierre Auger experiment both within the internal and external shock models. 
We have also studied the dependence of the maximum proton energies and the cooling 
breaks in the proton spectrum on the various parameters like Lorentz factor, energy of 
the GRB fireball, variability time (in case of internal shocks), ambient particle density 
(in case of external shocks). Our results are important to understand how the various 
observable parameters determine which mechanism e.g. pj interactions, synchrotron 
cooling of protons will dominate over one another inside these sources. 

PACS numbers:98.70.Rz, 98.70 Sa 
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1 Introduction 

The field of ultra-high energy cosmic rays (UHECR) enriched with the recent data from 
the Pierre Auger (PA) experiment [1] has emerged as a very exciting area of research. The 
highest energy cosmic ray event so far observed was of energy 300EeV, detected by Fly's Eye 
experiment [2] . Subsequently, many events of energies more than tens of EeV were detected 
by AGASA [3] and PA experiments. The origin of these extremely energetic cosmic ray 
events has been investigated since their detection. The current data from PA experiment 
shows anisotropy and the signature of GZK cut-off [4, 5]. Powerful cosmic ray accelerators 
e.g. gamma ray bursts (GRBs) and active galactic nuclei (AGN) within 200Mpc from us are 
likely sources of the events detected above the GZK cut-off energy [6] . A possible correlation 
of the AGN within 75Mpc with the cosmic ray events above 60 EeV has been identified by 
the PA collaboration [7]. Protons can be accelerated to above 100 EeV within GRBs both 
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in the internal shock [8, 9] and external shock scenarios [10, 11]. It was suggested by [8] that 
if all GRBs emit 5 x 10 50 erg in cosmic rays and the spectral index (p) of the accelerated 
proton spectrum (dN/dE oc E~ p ) inside them is p < 2 then they can explain the observed 
UHECR spectrum. The association of GRBs with UHECR has been further discussed in [9]. 
The energy production in cosmic rays was assumed to be same as the energy production in 
photons of energy between 20keV-2MeV which was used as 2 x 10 53 erg in this calculation 
and spectral index of proton spectrum p = 2. The local star formation rate was folded with 
the maximum distance from which GRBs can contribute to the observed proton spectrum at 
a particular energy and the proton spectrum from individual GRBs to calculate the diffuse 
proton spectrum. GRBs were identified as possible origin of ultra high energy cosmic rays. 
Within the external shock model it was shown that GRBs can explain the observed UHECR 
spectrum assuming all GRBs emit an average energy 3.3 x 10 53 erg in cosmic rays and proton 
spectral index p = 2.2. Although a steeper spectrum gave a better fit for the data below 
10 19 eV but due to the contributions of Galactic cosmic rays this fit can not be considered as 
better. Ultra-high energy (UHE) heavy nuclei can survive from photo-disintegration in the 
intense photon fields of GRBs in case of external shocks. It may also be possible in case of 
internal shocks if the shocks occur at relatively large radii [12]. 

A few low luminosity (LL) GRBs have been detected before e.g. GRB 060218/SN 2006aj 
and GRB 980425/SN 1998bw at redshifts of 0.033 and 0.0085 respectively. These GRBs 
show some common characteristics like long duration, low isotropic 7-ray energy and low 
luminosity. Their detections at very low redshifts within a comparatively short period of 
time imply that they are much more frequent than the canonical HL GRBs [13, 14, 15, 16]. 
The emission of UHE heavy nuclei and protons from high and low luminosity (HL and LL) 
GRBs have been studied by Murase et al. [17] considering internal and external shocks. HL 
and LL GRBs can emit UHE protons and heavy nuclei within the shock scenarios considered 
in their paper. Within the internal shock scenario LL GRBs are unlikely to produce 100 EeV 
protons but they can produce 100 EeV heavy nuclei. The values of the GRB parameters e.g. 
luminosity, Lorentz factor, variability time determine the maximum energy of the cosmic ray 
particles produced. 

The Galactic and extragalactic magnetic fields distort the trajectories of the charged particles 
and tracing the origin of these particles becomes very complicated. The magnetic fields are 
poorly known. For protons above 60 EeV the deflection in regular Galactic magnetic field is 
expected to be about few degrees [7]. For Fe nuclei of same energy the deflection would be 26 
times higher. The intergalactic magnetic field has been only measured at the centers of rich 
galaxy clusters. It can produce deflections from negligible [18] to very large values [19]. For 
coherence length of about 1 Mpc the intergalactic magnetic field may produce root-mean- 
square deflection of few degrees for protons of energy 60 EeV traversing through 100 Mpc 
[7]. As it is more difficult to track the heavy nuclei due to their large deflections, the protons 
of energy more than 60 EeV are only considered in this paper. Modeling of gamma ray 
bursts has remained as one of the most challenging problems in theoretical astrophysics for 
the last three decades. Although, lot of progress has been made in the past few years, as of 
now there is no unique model for GRBs which can explain all aspects and characteristics of 
GRB emissions in different wavelengths. The leading models for emission from these sources 
are the internal and the external shock models [20]. We calculate the proton fluxes from 
all GRBs within the internal shock scenario using distribution functions for their values of 
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luminosities, redshifts, Lorentz factors etc. The luminosity distribution functions of GRBs 
have been used from [16] and for redshift distribution we use star formation rate from [21]. 
The Lorentz factors, variability times and burst durations are represented by log normal 
distribution functions [22]. Considering the external shock model the time dependent and 
time averaged proton fluxes have been calculated. Our calculated fluxes have been compared 
with observed data on UHECR. 



2 UHECR Propagation and Spacial Distribution 

Cosmic rays produced from a point source are expected to reach us much later compared 
to the photons due to their deflections by the Galactic and intergalactic magnetic fields 
[23, 24]. If the cosmic ray protons are randomly scattered in magnetic field B of coherence 
length A then the typical deflection angle is 9 ~ (D/ \)°- 5 X/ R L , where D is the distance 
of the source, Rl = E p /eB is Larmour radius. The time delay of cosmic rays can be 
expressed as t(E) ~ D9 2 /4c. The dispersion in deflection angle A9 gives a dispersion in 
time delay of protons At(E p ). If the dispersion in deflection angle A9 << 9 then the 
dispersion in time delay is Ar(E) « r(E). In this case only a few HL GRB events are 
expected during At(E p ). The extremely energetic cosmic ray events from these GRBs will 
be spacially localized in a small region which is not consistent with the current data from 
PA experiment. The proton flux from a GRB in observer's frame is 

dN°\E° p b ) dN p (E p ) (1 + z)cos 2 9 a 
dEf dE p AttD 2 

x exp(-D/[cos9 a r n (E p ,z)]) (1) 

where Ef = E p /{1 + z) has b een used. We are assuming an average deflection of 9 a of the 
cosmic ray protons above 60 EeV. The comoving radial distance of the GRB is 

D = r±- , a (2) 

^o v /fi A + fi m (l + ^) 3 

with Vl\ = 0.73, Vl m = 0.27 and the Hubble constant H = 7lkmsec~ l Mpc^ 1 . Due to 
deflection in magnetic field the average distance traversed by the protons is D/cos9 a . The 
ultra-high energy (UHE) protons lose energy during their propogation due to photopion, 
photopair production. The mean free path of energy loss of an UHE proton with energy E p 
is [6, 25] 

r,(^)=(l + z)- 3 ^f^Mpc (3) 

where E Pi2 o = E p (eV)/10 20 eV . The proton flux from a GRB per unit energy and per unit 
time is 

dN°\E p ) dNp(Ep) (l + z)cos 2 9 a 

dE° b dtf ~ dEpAr(Ep) AnD 2 

x exp(-D/[cos9 a r\(Ep,z)]) (4) 

All the GRBs happening during At(E p ) contribute to the diffuse cosmic ray flux and the 
factor Ar(Ep) cancels out in the final expression of the diffuse flux. 
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3 The Different Time Scales, Cooling Break and Max- 
imum Energies of Protons 

We define the relevant time scales which can be used to calculate the break energy in the 
shock accelerated proton spectrum and the maximum energies of protons. The wind expansion 
time scale or dynamical time scale is defined as t' dyn = ^ where r is the radius of expansion 
and T is the bulk Lorentz factor. The synchrotron cooling time scale of protons in the 

comoving frame is t' = *L n , where the ratio of speed of protons to that of light is 

f3 p ~ 1, B' is the magnetic field in the comoving frame and a pT is the Thomson scattering 
cross section for protons. The cooling time scale due to J97 interactions is t' pi = Vcf r ^ E . From 
the condition t' dyn = t' syn we get the break energy in the proton spectrum due to synchrotron 
cooling E' pcsyn . The break energy due to py interactions E pcpi can be calculated from 
't'dyn = t'p-y The cooling break energy in the shock accelerated proton spectrum is close to 
the minimum of these two breaks 



E' pjc = min[E' p ^ syn , E' PjC J (5) 
The protons are shock accelerated in magnetic field by Fermi mechanism. Their acceleration 

o /- 771/ 

time scale is t' acc = -^f where C ~ 1. If we compare t' acc with t' dyn , t' syn and t' pj we get the 
cut-off energies in the proton spectrum E' pmaxdyn , E' pmaxsyn and E' p max pi respectively. The 
minimum of these three gives the maximum energy of the shock accelerated protons. 

Ep <m ax ^^^\^p,max,syni ^p,max,pyi ^p,max,dyn] (^) 



4 Proton Spectrum in the Internal Shock Model 

The internal shock model of prompt emission is widely discussed in the literature [26, 27, 28, 
29, 30, 31, 32, 33]. Its most successful aspect is its ability to model the complex temporal 
profiles of light curves during GRB prompt emission. The observed temporal behaviour 
reflects the temporal behaviour of the central engine. For long bursts the central engine 
is located deep inside a collapsing star, whose envelop may act as an additional agent to 
regulate the variability of the relativistic flow. Relativistically moving shells with Lorentz 
factors T ~ 100 — 1000 emerge from the core. The faster moving shell catches up the slower 
one. When the shells collide with each other their kinetic energy is transformed into internal 
energy of the fireball. Electrons and protons are shock accelerated in irregular magnetic fields 
when the shells collide. The total energy of the fireball is distributed among the electrons, 
protons and the magnetic field where e e , e p and e# denote the equipartition parameters 
respectively. The relativistic electron/proton spectra can be represented by broken power 
laws 

dNx(E') f E'- p E> x , m <E' X < E', c 

dE' x 1 ^r^ 1 e' x ,c<K [n 

in the case of slow cooling where x = e,p for electrons and protons respectively. E' 
denotes the minimum injection energy of electrons or protons. E' em = m e c 2 ^' p g{p)^-^- and 



4 



K, m = %rn p c 2 g(p) where g(p) = for p » 2 and g(p) ~ 1/6 for p = 2 [33]. % ~ 1. The 
cooling break energy is denoted by E' x c . We have derived its expression by comparing the 
dynamical expansion time scale (t' dyn ) with the cooling time scale (t' cool ). We use the general 
expression for f n (E p ) from [22] 



fn(Ep) = fo 



1 34^2-1 , Ep ^ a2 _! 

a 2 +l \E h ) 
1.34"!- 1 / E p \ Ql -l 



(8) 



p6 



Variability time, peak energy of the low energy (KeV-MeV) photon fluence in MeV and 
corresponding threshold energy of protons for interactions are denoted by t v , E ytPea k,MeV 
and Ept, respectively. a<i = (p + 2)/2 is the spectral index of the low energy photon spectrum 
above the peak energy E 1 ^ peak>MeV and below it ai = (p + l)/2 (slow cooling). 
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where r = T 2 ct v has been used. The magnetic field in the comoving frame inside HL and LL 
GRBs can be expressed as [34] 
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= 5 x 10 4 G 
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(10) 



where £ is the compression ratio, which is about 7 for strong shocks. Comparing t' dyn = ^ 



and t' 
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E the cooling break due to synchrotron emission can be derived 
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where (5' ~ 



From t'dyn = t'pj one § ets 
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The luminosity emitted in synchrotron radiation by electrons is 

Vete T 



1 + Y e 



Hot 



(13) 



L to t is the total luminosity of the GRB fireball. The radiation efficiency of electrons r] e = 
\{E' / E' e m ) 2 ~ p , 1] for slow, fast cooling respectively. The relative importance between inverse 
compton (IC) and synchrotron cooling of electrons is denoted by 



Y P = 



J e,IC 



J e,syn 



-1 + J I + 4r) e e e /e B 



(14) 



where L e ,/c and L etSyn (this is same as L 7 in our discussion) are the luminosities of radiations 
by IC and synchrotron mechanisms respectively. The proton flux in the source rest frame 
can be normalised as follows 



j-Ep 



E, 



dN p (E p ) 

'^KT dEp 



tp^tot 



(15) 



Ep^min = FE' pm is the minimum injection energy in the source rest frame. E tot is the 
total energy of the fireball in the same frame. The maximum energy of protons derived 
by comparing the acceleration and cooling (synchrotron and jry interaction) time scales are 
Ep tm ax,sym E pmaxpi respectively. 
Cc = 2j Sf = t'syn and C ~ 1 give 
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and t' acc = t' pi gives 



e = r 

- LJ p,max,p"/ J- 



1 + ai eB'cTt v / E' pb \ ai -i 



For, ai = 1.5 
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The maximum energy derived by comparing the acceleration and dynamical time scales is 



HL>LL eB'cTX 

p,max,dyn ^ 

_ fi v 1n l9 (^B~lL to t,5l) 1/2 

^p,max,dyn O X 11) eV 

1 2 



,x,dvn = U A 1U 

2.5 

,-LL , -,,,19 (^B~lLtotA7) 1/2 ^ n o , 

%,max,dyn ~ 1 X iU fT ( ! 9) 



One should also keep in mind that there are several limitations with the internal shock model. 
The thickness of the shell (^) may limit the maximum energy of the protons. Moreover, the 
relative kinetic energy between the colliding shells has to give the emission energy of GRBs. 
Hence, the radaition efficiency has to be high to attain the observed energy output during 
prompt emission of these sources. Also, the spectra observed by BATSE from some bright 
GRBs suggest at least within the same burst the peak photon energy distribution is narrow. 
This is hard to achieve within the internal shock model unless one invokes a strong bimodal 
distribution in the values of the bulk Lorentz factor. 



5 Proton Spectrum in the External Shock Model 

The expanding shells of GRBs push the external medium which results in external shocks. 
In [26] the authors proposed that thermalization of the spectrum of the GRB can be avoided 
by introducing a slight contamination of baryons within the original fireball. The baryon 
loading parameter is defined as 

Mc 2 V ; 

where M is the total baryon mass and E tot is the total energy of the photon-pair fireball. [35] 
(MLR) showed that for a certain range of values of T only a fraction of the initial energy 
is still in the form of radiation when the expanding fireball is optically thin. The rest of 
the energy gets converted into kinetic energy of the baryons. For high baryon loading the 
baryons can be accelerated such that almost all the initial fireball energy transforms into 
kinetic energy of the baryons. After an initial period of acceleration, the fireball expands 
freely. The energy in this phase is concentrated in a thin slab initially of roughly constant 
thickness 5r ~ r . After free expansion starts the shell thickness is given by [35] 

5r=\ Tr ° r s<r<r b 

[ r/r r > r b v ' 

where r s is the radius where free expansion starts and later the shell starts expanding linearly 
with r when r > r b . The baryons can radiate if their kinetic energy can be properly random- 
ized either by collisions with the external interstellar medium (ISM) or with slower/faster 
portions of the relativistic flow itself [36]. Finally, the fireball shell is decelerated by the ISM. 
During the initial interaction of the fireball shell with the external medium a reverse shock 
propagates into the fireball to stop it. A prompt fireball collects 1/T times its rest mass 
from the ISM by the time the fireball's radius attends the deceleration radius. The baryons 
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or cosmic rays can loose energy by photopion interactions p + 7 — > 7r°, ti + and synchrotron 
radiation inside the shell of thickness Sr. We derive the cooling break energy and the maxi- 
mum proton energy comparing the different time scale for the forward shock scenario. In the 
external shock model the expansion time scale at any time t is t' dyn = Tt. At the deceleration 
epoch it is where r dec is the deceleration radius and T dec is the Lorentz factor at that ra- 
dius. This radius can be expressed as r de ^ ~ (2.5 x 10 16 cm)(£ toti52 /n) 1/3 (1^^2.5) " 2/3 [20] and 
r dtc — (5-4 x 10 16 cm)(_E toti 5o/n) 1/3 (r (ieCi i)" 2/3 where E tot is the total energy of a blast wave. 
The comoving magnetic field at the deceleration epoch is B d ^ = cT dec (32imm p eB) 1 ^ 2 = 
39r deC 2.5(ne Bi _ 1 ) 1 / 2 G, B d( ^ = 1.23r (iec l (ne B _ 1 ) 1 / 2 G. n ~ 1 is particle number density of the 
ambient medium. 

For interstellar medium (ISM) the radius of wind expansion at any time t measured in the 
source rest frame is ([37]) 

r HL (t) = 4 x 10 16 (s to t i5 2/n) 1/4 th 1/4 cm 
r LL (t) = 1.25 x lQ 1% (E tot ^/nf'\ h l/i cm (22) 
The Lorentz factor decreases with time as follows 

r HL (t) = i9A(E tot , 52 / n y / \ 

T LL (t) = ll(E tot>50 /n) 1/8 t h 3/8 (23) 

where th is the time measured in the source rest frame in hours. We have used r = T 2 ct. 
The time dependent comoving magnetic field is 

B' HL {t) = 15GefE tot ^nVH h -*'* 
B' LL (t) = A2GefE tot ^'H h ^ (24) 

From t' dyn = t' syn , the break energy in the proton spectrum of GRBs due to synchrotron 
cooling of protons is 

E L p X syn = 6 x ltf 2 E tot ^n-^e B -H- h l/A eV (25) 
At the deceleration epoch this expression becomes 
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The time dependence of this cooling break energy is calculated below. The photon energy 
flux at the peak/break energy is proportional to t _1 , t~' p l 2 in the case of fast and slow cooling 
respectively. For p ~ 2 the two cases are almost similar. This is also true for the luminosity 
of the low energy photons if the threshold energy of photo-pion production is below the 
peak energy. If the threshold energy is above the peak energy then the luminosity would be 
proportional to ^-p/ 4 " 1 / 2 for both fast and slow cooling of electrons. The threshold energy 
condition for p7 interaction in the source rest frame gives -E-y^eaifc (t)E pb (t) = 0.3T 2 (t)GeV 2 . 
If it is slow cooling -E 7iPea fc(^) t^ 1 ^ 2 and in the case of fast cooling it is proportional to 
£-3/2 _ jf the threshold energy of p7 interaction is lower than the peak fluence energy in 
the afterglow photon spectrum then luminosity L(t) oc i _1 ,£ _P//2 in the case of fast and 
slow cooling respectively. For p ~ 2 we get L(t) oc t^ 1 in both the cases. Considering the 
external medium to be ISM fo L ' LL oc L(t) / \T A {t)tE lyPeak {t)) oc t°,t for slow and fast cooling 
respectively. The break energy due to p~y cooling in the source rest frame is proportional to 
t _1 / 4 ,t" 5 / 4 with cci ~ 1.5, for slow and fast cooling respectively. The time dependence of 
the p7 cooling depends on several factors like the photon flux, the external medium (ISM 
or wind) and whether electrons cool slow or fast. For different underlying assumptions one 
gets different time dependence of the p7 cooling ([12, 17]). The cooling break energy in the 
proton spectrum is E' = min[E' ,E' ]. We compare the cooling and acceleration 
time scales to derive the maximum energy of protons. 

rpHL 7 I r.21 dec,2.5 y 

^p^ax.syn 1 A iU (ri€ B 

E L p L max , syn = 4 x ^^^eV (28) 

After the deceleration epoch the time dependent maximum energy obtained by comparing 
synchrotron cooling and acceleration time scales is 

pHL _ r( +\ r 3em p c4 1 °- 5 _ f. , n 20 -1/4 pl/16 5/16 r 3/16 T/ / 9Q x 
&p,max,syn ~ 1 W ~ r>'HL(+\ ~ b X 1U € B E tot52 U t h eV (29) 

for f3' ~ 1. The maximum energy of protons derived by comparing the acceleration and pion 
production time scales is 



rpHL r eB' HL (t) T(t)t(l + a,) / ggV -ii v«i ^ 
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From t' dyn = t' acc we get 



F HL,LL _ eB HL ' LL (t)c 2 

p,max,dyn 

E"L X , d yn = 1-4 X lOi^^^LW^) 178 ^ 

E p,tax,dyn = 2-5 x 10 18 e B 1 /2 (£ ;3 oti5ori/ ^ ) i/8 ey (33) 
At the deceleration epoch 

rpHL,LL __ e -Pdec ffL 
p,max,dyn <^ ec 

C = 2 x 10 19 (r i e B ,_ 1 ) 1/2 (%^) 1/3 (r deCi2 . 5 )V3 e y 

E i,max,dyn = ^ X lO 18 ^,-!) 1 / 2 (^) ^OW) ^ ( 3 4) 

We calculate E Pymax as discussed in the previous section. The protons lose energy as long 
as the cooling time scale is much shorter than the wind expansion time scale. For e p ~ 1 
the mean Lorentz factor of the protons is roughly the Lorentz factor of the forward shock. 
The total number of protons increases as nr(t) 3 oc t 3 / 4 incase of ISM and as nr(t) 3 oc r(t) oc 
t 1 / 2 for wind medium. The total energy carried by protons remains constant, e p E tot = 
N p T(t) 2 m p c 2 oc t°. The proton spectrum has been normalised as given in eqn.(15). It is 
important to discuss the difficulty with the external shock model. The particles (e,p) need 
to be thrown back to the external shock region by some magnetic field for their reacceleration. 
The magnetic field in the ISM is known to be low. For acceleration of the protons to very high 
energy in the external shock model the ambient magnetic field has to be more than its typical 
value in the ISM. It is not physically impossible that the ambient medium close to GRBs 
is more magnetized than the average ISM. Alternatively, magnetic field amplification may 
lead to acceleration of particles to very high energy. For discussions on particle acceleration 
in this way in the non relativistic regime one may see [38, 39]. 



6 Internal Shocks:Average Proton Flux from GRBs 



We have used the broken power law luminosity distribution functions for HL GRBs [16]. 



L r L\(h 
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The break luminosity L b = 10 52 erg/sec for HL GRBs. The values of (3 2 are taken to be 
0.65, 2.3. The star formation rate is ([21]) 
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and p ~ 1. The diffuse proton flux can be calculated as follows 
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Figure 1: proton flux above 60EeV in eV 2 m~ 2 sec^ 1 sr^ 1 from HL GRBs using the internal 
shock model, z max = 0.02, 9 a = 4°, 10° (solid line, dashed line). We have considered fast 
cooling of electrons and p = 2.3 [33], e p = 0.8 and e e = 0.1. The data points are from recent 
results of PA experiment [43] 

Trf is the duration of a burst and t v is the variability time scale in the internal shock model. 
The isotropic energy emitted by low energy photons (E 1 = T d L 7 ) is related to the total 
energy of a GRB 

E i = Y+Y Etot ( 38 ) 

Amati-relation [40, 41] has been used to calculate the peak fluence energy in the low energy 
photon spectrum while using the internal shock model. We have averaged over all the relevant 
parameters (T, L 7 , Tj, t v ) as discussed in [22]. For the log normal distribution in Lorentz 
factor the peak is assumed to be at log{Y m ) = 2.6 with standard deviation 0.3 for HL GRBs. 
The durations of the bursts are assumed to follow log normal distribution with the mean 
at log{Td t m) = 1-5 with standard deviation 0.5. The mean of the log normal distribution 
in variability time is assumed to be at log{t v , m ) = —1.5 with standard deviation 0.3. In 
Fig.l. our calculated proton spectrum has been compared with PA data above 60EeV. If 
we increase the average deflection angle of the cosmic ray protons from 4° to 10° the proton 
flux does not change significantly. 

7 External Shocks:Proton Flux from GRBs 

In the external shock model we have calculated the proton flux from all GRBs using log 
normal distribution in total energy with the mean at log(E to t,m,erg) = 54 and standard 
deviation 2 [42]. The distribution in the values of the Lorentz factor at the deceleration 
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Figure 2: proton flux above 60EeV in eV 2 m~ 2 sec" 1 sr" 1 from HL GRBs using the external 
shock model at the deceleration epoch, z max = 0.02, 9 a = 4° as before also p = 2.3 [33]. 
e p = 0.5 and e# = 0.3. 

epoch is also assumed to be log normal with the mean 2.6 and standard deviation 0.3. The 
proton flux from all HL GRBs at the deceleration epoch is shown in Fig.2. We have used 
the particle number density in the ISM to be n = lcm~ 3 . After the deceleration epoch 
the maximum energy of protons decreases with time and the chances of getting cosmic ray 
protons above 60 EeV becomes less. From LL GRBs we expect even lower cosmic ray flux 
than that from HL ones. 

8 Maximum Energies of Protons: Internal vs External 
Shocks 

In Fig. 3. the three plots show how the maximum and cooling energies depend on the GRB 
parameters in the internal shock model. For T ~ 100 the minimum of the maximum energies 
of protons is due to jry cooling. As Lorentz factor increases the maximum energy of the 
protons is determined by the dynamical time scale. If the total energy of the low energy 
photons is varied synchotron cooling of protons remains unaffected as long as the luminosity 
of a GRB remains constant. For T = 300, t v = Isec and L 7 = 10 51 eg/sec the dynamical 
timescale determines the maximum energy of the protons in the entire range of E 1 considered 
in our work. The maximum energy determined by the dynamical time scale is independent 
of the variability time of GRBs as shown in the third plot of Fig.3. For higher values of 
variability times of GRBs the rate of dynamical expansion determines the minimum of the 
maximum energies of protons. Fig.4. shows the various energies (maximum and cooling) in 
the case of external shocks. Unless Lorentz factor is nearly 1000, pq cooling determines the 
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Figure 3: Internal Shock Model: First plot- total energy carried by low energy photons 
-^7,52 = 1, t v — lsec; second plot-r = 300, t v = lsec; third plot- T = 300, -E 7 ,52 = 1- E p ^ max , pi 
(solid line); E p>maXydyn (dotted line); E P)maXySyn (dashed line); E P ^ P1 (long dash dotted line); 
E PtC ,syn (small dash dotted line); e p = 0.8, e e = 0.1, L 1 = 10 51 erg/ sec 

maximum energies of protons. p~f cooling is more important in the case of external shocks 
than in internal shocks. We vary the total energy of the GRB, particle density of the medium 
and find jry cooling gives the minimum of the maximum energies of potons. The cooling 
break energies are not important in case of external shocks as they are very high. 

9 Conclusions 

The leading theoretical models often used to explain emissions from GRBs are the internal 
and the external shock models. Within these models cosmic rays can be accelerated to 
extremely high energies inside these sources. Cosmic ray proton fluxes above 60 EeV have 
been calculated from GRBs distributed upto a redshift of 0.02 assuming distributions in the 
values of the relevant GRB parameters. The recent data from the Pierre Auger experiment 
is compared with our calculated proton fluxes. Our results show that GRBs can emit cosmic 
ray flux comparable to the PA data above 60EeV both in the internal and external shock 
models. In the internal shock model p~f cooling and the dynamical expansion time scale of 
the GRBs are important in determining the maximum energy of the protons while in the 
case of the external shock model jry cooling determines the maximum energy of the shock 
accelerated protons. 
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Figure 4: External Shock Model: First plot- E to t,52 — 1, n — 1cm 3 ; second plot- Y = 300, 
n = lcm -3 ; third plot- T = 300, E to t,52 = 1; e.B = 0.5, e p = 0.3. Different line styles used as 
in Fig.3. 
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